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ff^ . Summary. We discuss the local behaviour of vector fields in the plane 



around a regular singular point, using recently introduced reduced normal torms, 
i.e. Poincare and Lie renormalized forms We give a complete clas- 

^ ' sification, and provide explicit formulas, using Poincare renormalized forms for 

non-degenerate cases, and Lie ones for certain degenerate cases. Both schemes 
are completely algorithmic, prove to be easy to implement, and only require to 
solve linear equations. 



1 Introduction 



O 

^ ; The theory and method of normal forms (NF) §, |, ||, ||, H ||, H , whose 

' origins go back to the work of Poincare at the end of XIX century, constitute 

a fundamental tool to study the behaviour of dynamical systems locally near a 
known solution. 

The normal form of a dynamical system x = f{x) (or equivalently of the 
' vector field X = f''{x)d/dx'^) with a fixed point in the origin depends on the 

I properties of the matrix A — (_D/)(0). As well known, when the eigenvalues are 

nonresonant the normal form is just given by the linear part of the system, and 
thus unique, while for resonant eigenvalues the normal form is in general not 
unique and can depend on infinitely many arbitrary constants. 

Needless to say, this richness of the normal forms unfolding reflects the rich- 
ness of possible different behaviours of nonlinear systems sharing the same linear 
part; however, it is well known that this is to some extent redundant: a given 
system with resonant linear part does not have a unique normal form, if not 
in exceptional cases (this is related to the freedom in choosing solutions to the 
homological equation at resonant orders k). Thus several authors have tried to 
devise ways to reduce this redundancy of the normal form classification, and on 
the other side to take advantage of the freedom mentioned above; in this respect 

*e-mail: g.gaeta@tiscali.it. 
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one should quote |, |, |, |T|, |T|, |o[ |3|, |o[ |5| as well as |4|, 
This problem was actually already mentioned by Dulac |26[|. 

One of these attempts, proposed in |3|], is based on considering higher 
order homological operators and the related homological equations (details on 
this approach will be given below); the resulting "further simplified" normal 
form has been called Poincare renormalized form (PRF). In a later con- 
tribution it was noted that the PRF approach does not make use of the 
Lie algebraic structure of vector fields in normal form, and that in certain cases 
one can obtain a stronger simplification by considering this structure^ The 
resulting "further simplified" normal form has been called Lie renormalized 
form (LRF). 

In this note I want to use PRFs and LRFs to analyze the behaviour of vector 
fields (dynamical systems) in the plane locally near regular singular points 
(equilibria). 

Together with general results, I will also give detailed computations up to 
some finite order (typically up to terms of order six in the x, y variables) with 
explicit identification of the transformations needed to take a system in PRF, 
including closed-form expression of the numerical coefficients. This will show 
that the required computations are actually easy to implement in practice. The 
analysis will be on the formal level only, which is standard in normal forms 
theory, and however useful for the analysis of the system in a way I will not 
discuss here (see e.g. |||,|3, 24, 53[ |5[ |2[ |3| for this matter). I will however show 



how one can explicitely determine the radius of convergence of the normalizing 
transformation at any finite order, see section 7.5. 

Some confusion arised in the literature about PRFs, in particular concerning 
a system considered by Bruno and Petrovich and the application of the PRF 
scheme to this; the appendix is devoted to clarifying this issue. 

We use frequently the abbreviations NF(s) for normal form(s), PRF(s) 
for Poincare renormalized form(s), and LRF(s) for Lie renormalized form(s). 
Equations are consecutively numbered in each section, and we omit the section 
number when referring to equations of the same section. 

The financial support of Fondazione CARIPLO per la ricerca scientifica 
during part of this work is gratefully acknowledged. 



2 Poincare normal and renormalized forms 

Let us start by introducing some general notation. We will consider vector fields 
in R^; we will use coordinates {x,y) in R^, corresponding to the basis (ei,e2), 
and denote a generic vector as f = (^^,^^). We will also write di = 9/9^*. 

We will denote by J- the set of polynomial vector functions, i.e. of polynomial 
functions / : R^ — > R^, having a zero in the origin; we denote by J^k C J- {k > 0) 
the set of polynomial vector functions homogeneous of degree fc + 1 in the i^. 

^ Doing sn is sn natural that the structure was inadvertently used in some examples con- 
sidered in teol, blj; this is the content of the errata to these papers; see also section 10 below. 
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We denote by W the Lie algebra of polynomial vector fields in equipped 
with the commutator operation. 

If we focus on the coordinate expression of vector fields, the role of the 
commutator is taken by the Lie-Poisson bracket {., .} defined as 

{f,g} := {fd,)9~{g^d,)f. (2.1) 

Indeed, writing Xf — pdi, Xg — g^di, we have 

[Xf,Xg] = . (2.2) 

The set equipped with the bracket {., .} is a Lie algebra. Notice that 
{., .} : J-j; J'm ^ ^k+m- 

We will also, with an abuse of notation, denote by Wfc the set of vector fields 
whose components are homogeneous of degree A; + 1 in the ^, and by Wk the 
homogeneous part of order fc + 1 of the vector field W . Obviously these are not 
intrinsically defined notions, but depend on the coordinates we use; thus if we 
consider a vector field W , when we change coordinates the Wk will also change 
(but near identity changes of coordinates — > C = + where ij^ S Tm, 

will preserve the Wk with k < m). 

To the linear part AS^ of a vector function f ^ T we associate the homological 
operator Cq — {A£^, .}. Notice that Cq : J^k ^ ^k- 

We can also define the homological operator in W rather than in JT, as 
follows. If the linear part of / is given by A^, we will denote the linear part (in 
the ^ coordinates) of Xf as Xa- To the linear part Xa of a vector field Xf we 
associate the homological operator Lq = [Xa, ■]', note that Cq : Wk Wk- 

We will equip Tk (and thus all oi V — Tq (B (B ....) with the Bargmann 
scalar product ||2^, |3^; this is defined as follows: 

With this choice^ of scalar product in V , the adjoint of is given by L\ = 
.}, where is the adjoint of A: A± = A*^. 
The operators Co and Cq play a crucial role in discussing the properties of / 
under Poincare transformations, i.e. under near-identity changes of coordinates 
in R2, given by C = C + ^HC), with h = hk G 

2.1 Standard Poincare normal forms 

It is well known that by a careful use of Poincare transformations, i.e. perform- 
ing them for k = 1,2,... successively and choosing the /ifc's as solution to the 
homological equations (see below), one can eliminate all terms in the range of 
Co. 

■^Using the standard scalar product |Q would differ, here and below, only in some coeffi- 
cients. 
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That is, one can pass to coordinates rj which reduce the coordinate expression 
of / to a form (the Poincare-Dulac normal form, or simply normal form) /, 
where /(??) = Arj + F{r]) and F e Ker(/:J). 

It is also well known that Ker(£j) belongs to (and coincides with for A 
semisimple) the set of resonant vectors, which are defined as follows. Consider 
a basis in such that A is in Jordan normal form, and let Ai,A2 be its 
eigenvalues (possibly equal). 

Then a resonant monomial vector ^^e^ of order \fj,\ is a vector v with 
components (we write the vector indices as lower ones for ease of notation) 
Vq = ^(i')?(2)<^a,/3 — x^^^y^^5a.i3, whcrc |//| = + /i2 > 1 and the are non- 
negative integers satisfying the resonance relation 

/iiAi + /i2A2 = \p . (2.4) 

The linear span of resonant monomial vectors is the space of resonant vectors, 
i.e. Ker(/:+)\[Ker(£+)n^o]. 

The homological equation for hk is given as follows: let / be the expression 
of / obtained after operating the previous Poincare transformations, and let tt^ 
be the projection operator nk : Ran(£o) H Tk', then 

Coihk) = TTkfk (2.5) 

is the required homological equation for hk', notice that the solution to this is 
uniquely defined up to elements of Ker(>Co). 

We refer e.g. to |, |,|4[|^, |3|||, || ||, |4| for further detail on standard 



normal forms and the normalizing transformation. 



2.2 Poincare renormalized forms 

In order to discuss PRFs, it is convenient to use Lie-Poincare - rather than 
Poincare - transformations. Let us first of all briefly discuss these, referring to 
e.g. §, H, H or ||, |l| for further detail. 

The function h : R^, or equivalently the vector field H = ft,'^' (x, y)dx + 

/i*^^-* {x, y)dy, generates a Lie-Poincare transformation given by the time-one flow 
of H. Thus under the Lie-Poincare transformation generated by the vector field 
H, the vector field W is transformed into 

W = e"We-"; (2.6) 

this can be computed up to any desired order by means of the classical Baker- 
Campbell-Haussdorff formula as 

W = J2-[[H,W]r, (2.7) 

where we have defined the iterated commutators as 

[[H,W]]° := W ; [[H,W]Y := [H , [[H,W]Y-'^ ] (s > 1) . (2.8) 
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If = /ife € ^fe) from the above we have, denoting by [a] the integer part of 
a and with 7i{f) := {h, /}, that 

[m/fc] ^ 

/m - 51 ^/^"(f'-n-sk) ■ (2.9) 

We can now define Poincare renormalized forms (PRFs) and describe their 
construction. 

We define the higher homological operators Ck as :— {fk, ■}', note that 
these make good sense only after fk has been stabihzed in the procedure, as 
discussed in [|o[|l). 

We define the spaces H^p'> Q T {p > 0) hy iJ^o) = T, and = n 

Kei{Cp) for p>0. This imphes that C and 

H^P'> = Ker(/:o)n...nKer(£„_i) = fX^eiiCs). (2.10) 

The restriction of jCp to H'-p^ will be denoted as ^Ap. With this definition, 
we have H^p+^^ = Ker(Xp). 

We also define the spaces F'^p'> <^ T {p > 0) us = T and F^p) = 

F(p-i) n Ker(X + ) for p>l. This implies that F^p+i) C F^p), and 

= fl' [Ran(A^,)]^ = fj' Ker(A^^) . (2.11) 

(the orthogonal complement must be understood in equipped with the scalar 
product). We also have F^^+i) = F^P\[RaTi{Mp) D F^p']. 

We can also define the projection operators TTfe : — > Ker(£fc), and lis = 
TTs-i o ... o TTo for s > 0, Ho the identity operator. Similarly, we define the 
projection operators Ps : T ^ Ran(y\/(s). Notice that with these H'^^^ = fipJ^, 
Mp = £p o ^p. 

The function / G (the associated vector field Xf ^ W) is said to be in 

FRF if fk e F^'^) (and then obviously f'' G F^''^ := F^'^) n J^fe). 

It can be shown that any function / G (any vector field W G W) can be 
taken to PRF by a sequence of suitably chosen Lie-Foincare transformations. 

Let us now briefly describe two possible schemes for constructing the se- 
quence of "suitably chosen" transformations; these were discussed in |3l| . 



(A). In the first case, denote by the term obtained after completing 
the procedure up to order k — 1. Then operate a series of transformations with 
generators^ h^i^\ h^kli, h[''~^\ with hp^ G H^'^^ HTp. These should be chosen 
as solutions to the higher order homological equations 

pji'^ - M. {hit) ; (2.12) 

•^The lower index keeps track of the subspace to which h belongs; the upper index keeps 
track of the number of transformations already operated. 
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in other words, 

h!^l - n,oX+on(/(^)) . (2.13) 

(B). Rather than putting /^"^ in F^"-* for s — 1,2, k, and doing this for 

(s) (s) 

all k = 1,2, we can invert the order of iterations, i.e. put ff, in Fj^ for 
all k > s, and do that for all s = 1,2, in this case for s = 1 we obtain the 
standard NF. 

Notice that the equations to solve, and the spaces to which the functions 
belong, are the same in the two cases; however, the form to which fk has been 

(s) 

taken by previous parts of the procedure when we deal with can be different. 

Due to non-unicity of PRF, these two procedure can indeed in principles give 
different PRFs, i.e. the arbitrary coefficient which appear in the general form of 
PRFs for a given system can take different values depending on the procedure 
we have followed. 

The reader is referred to [ p4| ^ for further detail concerning Poincare renor- 
malized forms and related matters, including the Hamiltonian version of the 
theory and the role of (linear) symmetries. 

We stress that the schemes mentioned here are "generic", i.e. do not take 
into account the Lie algebraic structure of the set Q of vector fields in normal 
form (with respect to a given linear part). This point will be considered in 
section 3 below, where a "^-adapted" procedure is discussed; this will also 
make transparent the relation between this approach and Broer's one []lT| 

3 Reduction of normal forms and Lie algebras 

The set of vector fields in (standard) normal form with respect to a given linear 
part has a natural Lie algebraic structure. It turns out that making use of this, 
one can obtain a better reduction of the normal form than using only the graded 
structure. 

3.1 The LRF approach 

Consider the set of vector fields in R" which are in normal form with respect 
to the given linear part A, i.e. the vector fields Y such that =0. It 

is obvious that these form a Lie algebra (the Lie operation being the standard 
commutator of vector fields); we denote this algebra by Qa- 

Let us recall a general characterization of vector fields in normal form rele- 
vant in this context |2j, |2^, Consider the linear vector field Xa; we say 
that the differentiable function (j> : R" ^ R is an invariant for Xa if XA{(f>) = 0. 

We will denote by X* (A) the set of invariants for Xa which are meromorphic 
(i.e. a quotient of algebraic functions) in the x coordinates; by T{A) C X*{A) 

would like to stress that the idea o f iis i pg Ck with the same role as Co was already 
contained in feol; at the time of writing ]30l Bl| I had not realized this, and did not give 
proper credit. 
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the set of algebraic invariants for Xa; and by Ik(A) C I{A) the set of algebraic 
invariants for Xa which are homogeneous of degree fc + 1 in the x variables. 

Let G = C{A) be the centralizer of A in the group GL(n); let its Lie algebra 
be spanned by matrices {Ki, Kj,} (we can always assume Ki = I, and that 
Ka = A for some a, provided A ^ 0; notice that d < n). The vector fields 
corresponding to these are given in the x coordinates by X^°'^ = {KaXydi. 

The most general vector field W in Q can be written as 

d 

where ^q(x) G X*{A); this implies that Q is contained in a finitely generated 
module over X*{A). 

As W must be algebraic in the x, and X^"^ are linear in x, functions ^a{x) G 
I* [A) with poles of degree d > 2 in a; = cannot appear in (10). In other 
words, only algebraic functions and functions with simple poles in the origin 
can appear in the actual normal form unfolding; hence Q is not, in general, the 
full G-generated module over I* {A) . 

In several cases, however, it happens that Q has a more convenient structure, 
i.e. the /Iq in (10) can actually be taken to be in I{A), and not just in X*{A). 
In this case we say that all the vector fields in Q are quasi-linear, or that we 
have a quasi-linear normal form. In particular, this is the case when A admits 
only one basic invariant (as in the cases we are interested in). 

If the normal form is quasilinear, we have Q n Vk+i = Ik{A) eg) G, and the 
analysis of the structure of Q results to be particularly simple. 

Call A'* the algebra spanned by vectors which are written as X = ^{x)X'^"'^ 
with ji Gl* (A) , and Xa the algebra spanned by vectors with ^ e X{A) (this is 
the module over X{A) generated by X^"^). In general we have Xi (B ■■■ (B Xd C 
Q C X^ © . . . © , and in the quasi- linear case we actually have Q = Xi® ■■■®Xd. 

It can happen that we are able to determine a sequence of subalgebras J-p C 

each of them being the union of Xa subalgebras, such that ^ Q and 

[g,Tp] = ^p+i ; (10) 

if this terminates in zero we say that Q has a quasi-nilpotent structure. Notice 
that the factor algebras Fp :— J-p/ J-p+i are in general not abelian. 

It results [ p2[ that Q can have a quasi-nilpotent structure only if G is nilpo- 
tent. The chain of subalgebras J-p (Z Q can then be read off the descending 
central series Gp of G. The subalgebras Fp introduced above are therefore mod- 
uli over X{A) generated by the abelian subalgebras 7p = Gp/Gp+i of G. 

If Q is quasi-nilpotent, we can first work with generators in Fi and simplify 
terms in Fi, then consider generators in r2 and simplify the corresponding terms 
being guaranteed that Fi terms are not changed, and so on. 

Notice that in this case we are - roughly speaking - just using the nilpotent 
structure of (the finite dimensional group) G, rather than the one of (the infinite 
dimensional algebra) Q. 
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Remark 6. It should be stressed that the reduced normal form obtain in 
this way is not necessarily a PRF in the sense of the definition discussed in 
section 2. We will use therefore the name Lie renormalized form (LRF) to 
emphasize that it is obtained using the Lie algebraic properties of the set of 
vectors in normal form and at the same time the main idea behind the PRF 
procedure. 

Needless to say, this approach is particularly convenient when the Tp are 
generated by a single element of G. 

The situation depicted above is met in applications: in particular, it applies 
to any nontrivial two-dimensional case. More generally, it always applies when 
there is only one basic invariant. 

A concrete application of this "^/-adapted procedure" will be given below 
when considering certain subcases, see subsections 6.3 and 7.4; in this case, 
indeed, the generic PRF procedure given in js^, |l] would produce an infinite 
PRF (as shown in subsections 6.1 and 6.2), while the ^/-adapted one will produce 
a finite LRF (as shown in subsection 6.4). In this case, it will turn out that the 
LRF is not a PRF. 

Remark 7. It should also be emphasized that this procedure can be seen 
as an implementation of Broer's idea on reduction of normal forms as filtration 
of Lie algebras; see also the work of Baider and coworkers. However, it is not 
exactly following their general procedure, based on descending central series (see 
the discussion in appendix B of [Q) but rather it is a blend of the Broer-Baider 
procedure and of the general PRF procedure. It should also be mentioned that 
the Brocr approach is extremely general, and produces unique normal forms; see 
[|[ for discussion and applications, and |45, 46| for a clarification of the 

issue. However, implementing it in practice can be of considerable difficulty, and 
as far as I know explicit computations based on this have been performed only 
for very simple systems. On the other hand, the LRF approach is not as general 
as Broer's approach (it requires the algebra Q to have a favourable structure) and 
will in general not produce unique reduced normal forms; but when applicable 
it only requires linear algebra computations and is easily implementable on a 
computer. 



4 Singular points of vector fields in the plane: 
the basic classification of linear parts. 

Let A = {Df){xo) be the linear part oi X = pdi at the equilibrium point xq; 
we can and will always shift coordinates in so that xq is in the origin. 

After reduction to Jordan normal form, and up to permutation of coordi- 
nates, the following cases are possible for A 7^ (all the constants /i, fii below 
are understood to be real and nonzero, and a — m + i/12): 
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In cases S1-S4 the matrix A is semisimple, in case Nl it has a seinisiinple part 
and a nilpotent one, in case N2 it is nilpotent, with zero semisimple part. Thus, 
case N2 (and the case A — 0) corresponds to a non- regular singular point ; it 
is known that in this case normal forms theory is not very effective and 
we will not deal with this; case N2 is studied from the point of view of PRFs 
in but the results that can be obtained are very poor (this singularity is 
better studied with different methods, see |Q and |3^). 



We will now briefly recall the results obtained by standard NF theory for 
each of the cases listed above. In several of them the PRFs are either trivial 
(i.e. coincide with the standard NF) or have been studied in [1,2]. 

In the generic case (SI) no resonance can be present (recall we assumed 
Hi 7^ 0), so the NF is linear; moreover the eigenvalues belong to a Poincare 
domain, so the normalizing transformation is convergent in some sufficiently 
small neighbourhood of the origin [|[ || . 

The case (S2), which is generic for hamiltonian system, has infinitely many 
resonances. The NF is written as W — [1 + a{\x\'^)]Ax + f3{\x\'^)Ex, where 
|xp — x\ + X2, a and (3 are arbitrary polynomial functions with zero constant 
part, and E is the identity matrix; in the hamiltonian case we obviously have 
/3 = 0. The further reduction of this NF has been considered by Siegel and Moser 



in the hamiltonian case (see |28 for higher dimensions), while the generic 



case /3 ^ has been studied via PRFs in ||30|, 31 (see section 10 below). Note 
that for /3 = the NF satisfies "condition A" and thus, provided the linear part 
satisfies the arithmetic "condition w" |lj, Q , the normalizing transformation is 
convergent on the basis of Markhashov-Bruno-Walcher-Cicogna theory [|8[ ^ 
^ I39I , [ssl ; no convergence result is available in the generic case. 

In case S3 the eigenvalues cannot belong to a Poincare domain, and it is 
easy to see that the NF will depend on two infinite sequences of real constants, 

k=l '^kX , N 

I, fc \^-^) 

bkX^y . 

The PRF in this case will be studied in detail in sections 6 and 7 (see also the 
appendix) . 

In the case S4 we should distinguish several subcases according to two cri- 
teria: first, if /ii//i2 is rational or irrational, and second according to the sign 

of ^1^2- 

For fiifi2 > the eigenvalues belong to a Poincare domain, and the trans- 
formation to NF is guaranteed to be convergent on the basis of the Poincare 
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criterion; if is irrational, the NF is linear, otherways it can include reso- 

nant nonlinear terms (see subsection 8.1 for limitations on these). 

For /Lti/Lt2 < the eigenvalues are not in a Poincare domain and we are not 
guaranteed of the convergence of the normalizing transformation on the basis 
of the Poincare criterion. If fXi/ fi2 is irrational, the NF is linear, and no further 
normalization is needed; moreover, convergence can be guaranteed on the basis 
of Pliss theorem Hi. If jii/ U2 = p/q & Q, there can be resonances; in this case 
Sternberg theorem^, |8[ |9|, ^ guarantees that the NF is smoothly (but in 
general, not analytically) equivalent to the original system. 

The PRF for this case has not been studied so far, and will be considered 
below in section 8. 

In case Nl the eigenvalues are in a Poincare domain, as /i ^ 0; moreover 
there are no resonances, and thus we have a linear NF, with a convergent nor- 
malizing transformation. 

Finally, in the nonregular case N2 the standard normal form is given by 



. k 



EcxD k , 1, k+1 (4-2) 

As already mentioned, the PRF for this case was studied in detail in |^. We 
briefly recall the (poor) results concerning this in section 10. 

It follows from the above summary of known results that we need to study 
PRFs only in the cases S3 and S4 (and, as mentioned in remark 1, to correct 
a formula in case S2). Actually, as discussed in the previous section, formal 
computations in one of these cases can be mapped to other cases as well. We 
will thus study the case S3 in full detail. 



5 The S3 case: standard normal forms 

Let us consider a linear part given by 

i.e. corresponding to the vector field ydy. We note immediately that here 
A = A"*", so that the homological operator associated to /o = Ax will satisfy 
Cq = Cq (recall we arc using the Bargmann scalar product). It is easy to see 
that the kernel of Cq is spanned by the arrays of vector fields (with A; > 0) 

Xk-.^x'^+^d, eWk and Yk:=x''ydy eWk (5.2) 

(with this notation the linear part considered here is given by Iq)- These vector 
fields satisfy the commutation relations 

[Xk, Xm] ^{m-k) Xk+m , [Xfe, Yrr,] = mYm , [Ife, Y^] = . (5.3) 
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Wc denote by Q the (infinite dimensional) Lie algebra spanned by the X^'s and 
the Yfe's; and by X the algebra spanned by the Xj^s, by y the algebra spanned 
by the Y^s. Obviously Q = X (B y- Note that y is an abelian ideal, actually 
the maximal abelian ideal, in G. The (standard) normal form corresponding to 
the linear part considered in this section will thus be given by a vector field 

oo 

W = Yo + ^(afcXfc + bkVk) (5.4) 

k=l 

depending on the two infinite sequences of real constants am,bm- This is pre- 
cisely the structure considered in section 3. 

Remark 9. The vector fields Z- := dx, Zq := xdx and Z+ := x^dx act 
in each of X and y, as respectively a lowering operator, a counting one, and 
a raising one: that is, = (n + l)X„_i, [Z_,Yn] = nYn-i; = 

(n+ l)X„+i, [Z+,Yr^] = nYn+i; [Zo,X„] = {n+l)Xn, [ZQ,Yn] = nY^. 



6 The S3 case: Poincare renormalized forms 

We want now to consider the PRF corresponding to the linear part given by 
(5.1). In the spirit of PRF, we should act on the NF (5.4) with Lie-Poincare 
transformations generated by homogeneous functions /im e Ker(£o)nym. These 
will correspond to the action of vector fields of the form iJ^ = aX^n + /S^m- 

We will first, in this section, consider the spaces defined in the PRF pro- 
cedure, and thus obtain the general form of the PRF in this case (in the next 
section, we will perform explicit computations up to order N = 5). We recall 
that Ker(£o) corresponds to the sum of the algebras X and y, and that here 
Ker(£o) =Ker(/:+). 

Wc have then to consider Ci ; this depends on the coefficients of the quadratic 
part Wi of the vector field W, which we write as 

Wi=aiXi+biYi . (6.1) 

The cases to be considered are 



ai ^ 0, 5i = 
ai = 0, 5i ^ 



(a) 



ai = , 6i = . (d) 

We refer to cases (a),(b),(c) as nondegenerate (although properly speaking 
only (c) is such), and to (d) as the degenerate (properly speaking, completely 
degenerate) case. 

6.1 The nondegenerate cases 

In case (a) we have Wi = aiXi; we notice that 

[Xi,Xk] = ik-l)Xk+i , [Xi,Yk] = kYk+i (6.3) 
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and therefore for Aii the restriction of Ci to Ker(£o) we have that Ker(A^i) 
reduces to the hnear span of {Fq)-''^!}, i-e- of Wq and Wi, so no further nor- 
maUzation employing Li, L^, ... is possible. We also have that the range of A^i 
(the most relevant space for our discussion) is the whole linear span of the {X/;} 
(with k>2) and of the {Yfe} (with k>2). 

Remark 10. We stress that for the sake of the present computations 
(which aim at identifying linear subspaccs) we can as well assume ai = 1; the 
same remark would apply to other cases. Such a trivial remark will be of use 
later on in section 8. 

In case (b) we have W\ = h{Y\. We notice that 

[yi,Xfe] = -n+i , [yi,yfe]=o, (6.4) 

and therefore Kcr(Ali) = y. On the other hand, Ran(A^i) also is given by y, 
and Ker(7Wf) = X. 

In this case we also have to consider higher order parts of W\ the first step 
of the PRF procedure can eliminate all terms in Ran(A^i) and thus we will only 
consider terms in Ker(Ali'). Let p be the first integer for which ^ 0, and let 
Wp = ttpXp (all the Yk parts with k > 2 can be eliminated, as just recalled). 
Now A4p is the restriction of Cp to Ker(7Wi) = Ker(£o) H Ker(£i): indeed the 
£m with 1 < m < p are zero and put no restriction. We have 

[Xp,Yk] = kYk+p (6.5) 

and thus Key{Aip) = {0}: no further normalization is possible. 

In case (c) the situation is quite similar to the one met in case (a): we have 
indeed Wi = aiXi + biYi with nonzero constants ai,6i; we have immediately 
that 

[Wi,aXk + l3Yk] = aia{k - l)Xk+i + (fcai/3 - bia)Yk+i (6.6) 

This shows that Kcr(Ali) is just given by {lo)Oi-^i + ^i^i}) i-e. by the linear 
span of Wo and Wi : so again no further normalization using operators £2, 'Cs, ... 
is possible. As for Ran(Ali), this is the linear span of {Xk\ with fc > 2, and of 
{Yk] with k>2. 

We summarize the results of this discussion as follows, with W the vector 
field after the whole PRF procedure and omitting the degenerate case (d) . The 
hat on constants a/j will indicate that coefficients are not the same as those of 
the initial NF (5.4). 

_ (Yo + axXr+a2X2 (a) 
W=IYo + hYi + Er=2 akXk (b) (6.7) 
[ Yo + aiXi + biYi + 02X2 (c) 

We anticipate that a different reduction scheme (see section 3) can give a 
finite dimensional NF in case (b); we discuss this later on in subsection 6.3. 
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6.2 The degenerate case 

The discussion of the degenerate case (d) would require to consider the first q 
such that Ug + bg ^ 0, and then repeat the considerations presented above in 
cases (a),(b),(c), obviously with the role of ai, 61 taken by aq,bq. This is done 
in the following lines. 

We denote by /U > 1 the first k such that Uk ^ 0, and by > 1 the first k 
such that bk 7^ 0; at least one of these has to exist and be finite, or the system 
would already be linear and thus trivial. 

We will have to consider three cases 

H <v (da) 

11 >v (db) (6.8) 
11 = u (dc) 



In case (da) the NF will be given by 

u—l 00 

W = + ^afcX^ + Y.^akXk + buYk) . (6.9) 

fc=/x k=v 

We write again a, Hk E Ker(£o) H Wk as Hk = OkXk + fhYk, and we have 

£^{H^) := [W^,Hk] = a^[X^,akXk + PkYk] = 
= a^a{k - n)Xf^+k + a^f]kkY^+k ■ 

Thus it suffices to operate successively transformations generated by Hk 
(with fc = 1, 2, ...) and choose at each step 

ak = > = , (6-11) 

[k - iijttfj, kttfj, 

where a^+kib^+k denote the coefficient of X^j^ki^mu+k in W, i.e. after the 
action of previous transformations. Notice that in this way we can eliminate all 
terms except the ^2,^ one {k = fi). Thus, the PRF in case (da) results to be 

W = Yo + a^,X^u + -qX-i^, (6.12) 

where is the same as in the NF and is a real number. 
In case (db) the NF is 

1 00 

k—v k—v 

Now we have for C^{Hk) that 

U{Hk) = b,.[Y^,akXk + l3kYk] = -ub^akY^+k (6.14) 
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and therefore we ean eliminate all the K^+fe terms simply by choosing, with the 
same notation as before, 

ak = ^ ; (6.15) 

we cannot eliminate any of the Xk terms. 
Thus, the PRF in case (db) is 

oo 

W = Yo + KY, + ^SfeXfe . (6.16) 

Similarly to what happens for the nondegenerate case (b), the reduction 
scheme based on LRF, discussed below, gives better results in this case. 
In case (do) we have ij. = f, the NF is 

oo 

W = Yo + Y,iakXk + bkYk) . (6.17) 

k=iJ. 

In this case 

jCi^{Hk) = [a^X^ + bf,Y^,akXk+ PkYk] = 

= a^afe(fc - ^)X^+fc + {kan(jk- lJ.h^ak)Y^+k ■ 

Thus for k Id it suffices to choose 

ak = = {k - ^i)a^h^+k + lib^^an+k 

{k - iJ,)af^ ' fc(fc - fi) 

to eliminate both the X^^^k and the Y^+fc terms. For k = fi, wc choose = 
and (3k = b^+k/ {ka^) and eliminate the Y2^ term. Thus, the PRF in case (dc) 
is 

W = Yo + a^X^ + b^Y^ + rjX2^, . (6.20) 

6.3 LRF reduction scheme for cases (b) and (db) 

In the previous computations, we have followed the general PRF scheme for 
further normalizing the standard NF (5.4); this gave an infinite PRF in case 
(b) and in the corresponding degenerate case (db). However, as discussed in 
subsection 3.3, one can take advantage of the specific Lie algebraic structure 
oi Qa ^ X ®_» y (this recalls that X is acting on y by inner automorphisms) 
to obtain a more drastical reduction: indeed, in this case the Lie renormalized 
form will be a finite one, as we now discuss. We use the same notation as in 
discussing the case (db) above. 

We first operate a sequence of normalizations with generators /i^"-* = akXk, 
which we choose so as to eliminate higher order Xk terms, i.e. Xk for k > jjl 
(as we know, this is not possible for k — 2/i). Notice this will change not only 
the (coefficients of the) Xk terms, but the (coefficients of the) Yk terms as well; 



14 



however, no terms of degree k < v will be produced. In this way, we arrive at 
a form of the type (as usual the tilde indicates that the cocfhcients are not the 
same as the initial ones, but not yet final) 

oo 

Once this has been done, we pass to consider a second sequence of normal- 
izations with generators h^^'^ — PkYk- As y is an ideal in Qa, in this way only 
Yfc terms are generated, i.e. the Xk terms are unaffected. On the other side, y 
is abelian, and so only the Xf^ and X2^ are actually active in these transforma- 
tions: that is, we can only eliminate terms Y^+i and higher (it is clear by the 
commutation relations that these can always be eliminated). 

In this way, we arrive at the LRF: this is a NF depending on — v + 3) 
constants (this agrees with the number of constants predicted by Bruno in sect. 
III.2.3 of Q; see also (H), of the form 

W ^ Yo + a^X^ + a2^X2^ + J^hYk . (6.22) 

k=i' 

It is also clear by this discussion that bk = b^, 0,2^ — 0-2^- 

Remark 11. It should be stressed that the LRF (22) is not a PRF. Indeed, 
in this case the spaces v < k < ^ reduce, as seen 

in subsection 6.2, to multiples of Xk- Here we have therefore Wk ^ i^^'''* for 
V < k < fi, and thus (see section 2) the LRF cannot be a PRF. 

7 The S3 case: explicit computations 

As stressed in the PRF procedure is completely constructive; indeed, the 

PRF procedure gives an algorithm (which is easy to implement on a computer, 
as shown by the formulas reported in this section) to determine the coefficients 
ak,Pk of transformations needed to take the system (5.4) into its PRF. 

In this section I follow these computations (not needed if one is only inter- 
ested in the most general PRF form) for the linear part A given by (5.1), up to 
order six in x and y - i.e. put W in PRF up to terms W5 - for the nondegen- 
erate cases.| To avoid trivial steps, I will always assume that a first Poincare 
normalization has already been performed, taking the system into its standard 
normal form /(i' {W = W^^'^). 

In order to display the rather long explicit formula we obtain, it will be 
convenient to denote by Wk the part of the (coordinate expression of the) vector 
field W homogeneous of degree A: -I- 1 in the coordinate we are using; note they 
will change with the changes of coordinates. 

^Explicit formulae for a specific degenerate case which raised some confusion in the Htera- 
ture are given in the appendix, up to order ten in x and y. 
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Recall the explicit expressions given below are computed using the Baker- 
Campbell- Haussdorff formula (1.7), (1.9); i.e. by considering Lie-Poincare trans- 
formations, and not simply Poincare ones. 



7.1 Case (a) 

We first operate a transformation with hi £ Ker(£o) H J-i, i.e. with Hi = 
a.\X\ -\- PiYi G Wi n Ker(£o); after this, the quadratic part Wi of the vector 
field is unchanged, W^^^ = wl^\ while the cubic one is given by 

^ = 02X2 + (62 - aiPi)Y2 . (7.1) 

We know from our previous general discussion that - as indeed obvious from the 
above formula - the first component of this cannot be eliminated; to eliminate 
the second one, we have to choose 

Pi = 62/ai ; (7.2) 

we can choose ai as we like, say 

ai = (7.3) 

for simplicity. This choice of ai, f3i fixes the PRF after the first renormalization, 

(2) 

i.e. /j; . We do not give the explicit formulae for the sake of brevity. 

Let us now operate a transformation with /i2 ^ J-2 r] Ker{Co); the terms 
/o, /i, /2 are unaffected. Using the explicit formulae for fj^ , we have that W3 ' 
is changed into 

WP = (03 - aiaa) X3 + (63 - 2oi/32 - 0262/01) F3 • (7.4) 

We know from our previous discussion that we should be able to eliminate both 
components of this vector; this can indeed be obtained by choosing 



"3 „ {aiba - 0262) , , 

"2 = — , P2 = TT^ (7.5) 

ai ia\ 



This choice of a2,P2 fixes the PRF after the second renormalization, i.e. the 
fj^^ . Again we do not give the explicit formulae. 

Let us now operate a transformation with /13 E J^^ n Ker(£o): the terms 
/o, ...,/3 are miaffected. Using the explicit formulae for f^^ and fj^'^\ we have 
that VF4 is changed into 

W]^) = {ai - 2aia3) ^4 + 

+ ([o2^2-oia2?>3+oi(-a3&2+oi(&4-3ai/?3))] /(a?)) I4 

(7.6) 

Again we know apriori that this can be eliminated, and indeed the above formula 
shows that this is the case if we choose 

as = a4/(2ai) 

03 = [02^2 - 010362 - 010263 + o?64]/(3of ) ^ ■ 
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This choice of 0:3 , (3^ fixes ff, ' . 

Let us now operate a transformation with /14 G n Ker(£o); the terms 
fo,---, fi are unaffected. Using the expUcit formulae for fj^\ fj,^^ and fjf'\ we 
have that Wg^^ is changed into 

W^^^ = [{al - 0204 + 2ai(a5 - 3aia4)) / (2ai)] X5 

+ [(a? (-0462 + 0363 - a2&b + 0165 - 4alPi) (7.8) 
-alb2 + aialba) / (a?)] Y5 

which goes to zero if we choose 

a4 = (al — 0204 + 2aia5) / (Oaf) 

(3i = - {alb2 + 0^0462 - ciialbs - afasba + 0^0264 - afb^) / (4af ) . 

(7.9) 

Clearly, the computation could be performed up to any desired order, com- 
patibly with the computational power at our disposal, producing more and more 
complex but still completely explicit formulae; we will stop at this order. 

7.2 Case (b) 

Let us now consider the (slightly more complex) case (b). This was subject to 
some controversy (see remark 2 and the appendix), so that we will discuss it 
in full detail, following the transformation of the coordinate expression of the 
vector field W step by step. 

With a first transformation generated by Hi = a\Xi + (}\Yi we have that 
W2 is given by 02X2 + (62 + ai6i)i2; requiring the coefficient of I2 to vanish, 
we get 

ai = , (7.10) 

and for the sake of simplicity we will take /3i = 0. In this way we get (here and 
in the following 0(6) denote terms in We and higher) 

W?(2) = + biYi + a-2. X2 + 

+ [as - a2&2/&i] + [63 - bl/bi] + 

+ [a4 - 20362/61 + a2bl/bl] X4 + 

+ [64 + 253/52 - 35253/61] I4+ (7.11) 
+ [as - 80462/51 + 3a35i/6f - 026^/6?] X5+ 
+ [65 - 354/5? + 65253/52 - 45254/5i] 
+ 0(6) . 

We will now operate a transformation generated by H2 = 02^2 +/32i2- This 
leaves lower order terms unaffected, while W3 reads after this 

^ = (as - (x4a262)/5i)X3 + {{-{{xSjhl) /b,)) + xSjb:, + xSMa2))Y^ (7.12) 
Requiring the vanishing of the coefficient of the Y^ term we get 

a2 = {bj - b^bs) I {b\) ; (7.13) 
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we will set again /32 = for the sake of simplicity. With these, we get 

= Yo + biYi + 02X2 + [a3-a2b2/h]Xs + 
+ [a4 - 20362/^1 + a2bybl] X4 + 
+ [2bl/bl - 36263/61 + 64] Y4 + 
+ [05 - 30462/61 + 3036^/6? - a26i/6?+ 

+ {a,b, - a2b2){bl - bib,)/bl] X5 + 
+ [96^/6? + 96163/6? - 36i/6i - 46264/61 + 65] Y5 + 0(6) . 

(7.14) 

Let us now consider a transformation with generator H3 = 0:3X3 + 
Now lower order terms are unaffected, while we get that the coefficient of the 
Y4 term is changed to 

i2bl)/{bl) - (36263)/(6i) + 64 + 6ia3 . (7.15) 

Requiring this to vanish, we get 

as = - (26^ - 3616263 + 6^64) / (6?) . (7.16) 

We will, as by now usual, set /?3 = 0; with these we obtain 

= Yo + biYr + 02X2 + [03-0262/611X3 + 

+ [04 - 20362/61 + 026^/6?] X4 + 

+ [05 - 30462/61 + 4036^/6? - 0363/61 + , . 

-2026263/6? + ^264/61] X5 + ^'--"^ 

+ [-961/(26?) + 96163/6? + 

-36i/(26i)- 46264/61+ 65)] + 0(6). 

We now operate with H4 = 04X4 + (34Y4; the coefficient of Y5 results to be 

-(96|)/(26?) + (96i63)/(6?) - (36i)/(26i) - (46264)/(6i) + 65 + 61^4 . (7.18) 

Requiring this to vanish, we get 

"4 = (96| - I8616I63 + 36?6i + 86?6264 - 26?65) / (26f) ; (7.19) 

we also set ^4 = 0. We now have 

M?(5) = Yo + biYi + 02 X2 + [03 - 0262/61] Xs + 
+ [04 - 20362/61 + 026i/6?] X4 + 

+ [05 - 30462/61 + 4o36i/6?+ ^ ■ ' 

-0363/61 - 2026263/6? + 0264/61] X5 + 0(6) . 

Again we will stop at this order; the result of this explicit computation fits in 
the general result obtained in the previous section. 
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7.3 Case (c) 

Wc could analyze the other case (c) and produce explict formulas proceeding 
in the same way as in the previously considered cases (a) and (b); however, the 
procedure is by now clear and for the sake of brevity we will just give the final 
formulae. 

The coefficients a are chosen as 



ai = 
012 = as/ai 

0:3 = a4/(2ai) 



(7.21) 



"4 = (ai - 01204 + 2aia5) / (6a^) ; 
the coefficients /3 are chosen as 

01 = &2/ai 

02 = (0361 -0262 + ai&3) /(2a?) 

03 = — (2a2a36i — aia4bi — 2a^b2 + 2(11(1:^1)2 + 2aia2&3 — 20^64) / (6af) 

04 = (Sa^asbi — 010361 — 201020461 + 010561 — 80262 — 80^0462+ 

+8010^63 + 80^0363 - 80^0264 + 8of 65) / (12of ) . 

(7.22) 

It this way, wc arrive at a PRF given by 

= Yo + oiXi + 6iyi + 02X2 + 0(6) . (7.28) 

This corresponds to the result obtained by our general discussion above (more- 
over the coefficient of the X2 term is unchanged). 



7.4 The LRF scheme for ceise (b) 

As mentioned in subsection 6.8, in case (b) the alternative scheme adapted to 
the structure of G described there permits to obtain a finite NF (the LRF) and 
is thus to be preferred to the general one. Here we briefly discuss the explicit 
computation to be performed according to this. We deal with the nondegenerate 
(properly speaking, not completely degenerate) case, which means fj, = 2, p = 1; 
see subsection 6.3. 

With a transformation hi = aiXi, the W3 term reads 

W3 = [03 + oaai] X3 + [63 + 262^1 + 61a?] F3 . (7.24) 

We disregard the Y3 term and choose ai so to eliminate the X3 term, i.e. 
ai = —03/02. After computing the effect of this on higher order terms, we 
could perform a transformation with generator /12 = a2X2. However, we know 
that there will be no way to eliminate the X4 term, so we set a2 = 0. We 
perform a transformation with generator hs = 03X3. With this, the W5 term 
reads 

W5 = [2o|/o| - 803O4/02 + 05 - 02a3]X5 + 

+ [ajbi/a^ - 4a|62/a^ + 60^63/01 - 40364/02+ (7.25) 

+ 65 - 203610:3/02 + 262Q!3]l5 . 
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Again we only aim at eliminating the X5 term, and thus we choose 

2a| - 8020304 + al 05 

as = -3 . (7.26) 

"2 

We will be satisfied with this order of normalization for the X). terms, and take 
now care of the Y/. ones. 

We operate a transformation with generator hi = with this we have 

that _ 

W3 = [al5i/a2- 20362/02 + 63 -a2/3i] I3 • (7.27) 

We eliminate this by choosing 

alh - 2020362 + alb3 

Pi = -3 • y'-^^) 

"2 

We compute the effect on higher order term, and then consider a transformation 
with generator /12 — /32i^2; with these, we have 

W4 = [04 - a§/o2] X4, + 

+ (1/oi) [api + 3020^62 - 80203(0461 + 0263)+ (7.29) 
+0^(0561 + 02(64 - 202/32)) ] Y4 . 

We want to eliminate the Y4 term, and thus we choose 

f32 = —-4 (o|6i - 802030461 + O2O561 + 802O362 — 802O363 + O264) (7.80) 
202 

Again we take into account the effect of this on higher order terms, and pass to 
consider a transformation with generator /i4 = /34I4; we get 

W5 = [—20161/02 + 5030461/02 — 2030561/02 — 20162/02^4030462/01+ 
+20562/02 + 7o§63/a| - 0463/02 - 40364/02 + 65 - 802/^3] Y^ 

(7.31) 

which can be eliminated by choosing 

133= - (1/(802)) (20361 - 502030461 + 203030561 + 2020362+ 

+403030462 — 2020562 — 7030363 + 030463 + 4030364 — 0365) . 

(7.32) 

Summarizing, and having taken into account all higher order effects (up to 
order six), we have reached the LRF 

W = yo + 6iYi + 02X2 + [62- (0361/02)1^2 + [a4-(ai/a2)]X4+ 0(6) (7.88) 
We will be satisfied with this order of normalization. 
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7.5 Changes of coordinates 

In this section wc have given completely explicit formulas for the generators of 
Lie-Poincare transformations and for the PRF which can be obtained in this 
way for case S3 (these are also applied to other cases, as discussed in section 3 
and also in sections 9 and 10). It should be mentioned that in this simple case, 
one can describe exactly the change of coordinates generated by the vector field 
Hk = akXk + PkYk- 

The evolution under the vector field Hk {k > 1) is described by 

J=a,x'=+i , f^=Pkx''y (7.34) 

with initial datum x{0) = xo,y{0) = yo. The first of these is solved by elemen- 
tary methods to give 

^(^) = n ? kwk ■ (7-35) 

(1 — akksx^)'-'^ 

Using this expression for x{s), the second of (7.34) is rewritten 

^ = ,^ k, ds , (7.36) 

y (1-afcfcsxg) 

which gives 

y^'^ = ^° (1 - a.fc4)/5^/(-^) ^'-''^ 

We are interested in the mapping {x,y) = {xo,yo) (2;(1),2/(1)) := {x,y), 
and from the above we have that 

_ X ^ 1 

^ " {1 - akkx'^y/'' ' y = y {1 - akkx'^ f^/C'^^) ^^"^^^ 
with the inverse change of coordinates given by 

a^x ^ ' ^ 



1 + QffeX \ 1 + QffeX / 

These allow to obtain explicitely the changes of coordinates performed in 
passing from NFs to PRFs in case S3, and can be mapped to consider the other 
cases as well. The explicit formulas, however, would contain rational power and 
be quite involved. 

In order to make contact with the explicit formulas given above, notice that 

if we are acting with Hk , wc arc actually considering the change of coordinates 
from x'^ to a;''^+^). This map is defined only for x^ < (afcfc)"^; this allows 
therefore to explicitely compute the domain of analyticity of the change of co- 
ordinates. 
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8 The S4 case: standard normal forms 



We consider now the case S4, i.e. the hnear part of our vector field is now given 

with X ^ fi, X ^ and ^ 7^ 0. 

As remarked in section 4, if A/i > (the eigenvalues have the same sign), 
we are in a Poincare domain, so the convergence of the transformation to NF 
is guaranteed; on the other hand, if A/i < (i.e. we have an hyperbolic saddle 
point in the origin), we are not in a Poincare domain. However, the Chen- 
Sternberg theorem |^, |[ ^ guarantees the system is C°° conjugated to 
its normal form; as for the analytic conjugacy in this case, this is guaranteed if 
|A//i| is irrational, due to Pliss' theorem |Q. 

We also noticed that if X/ fi is irrational, there are no resonances, i.e. the NF 
is linear; in this case we do not need (nor it makes sense) to consider PRFs. Let 
us thus focus on the rational case. We will assume |A//i| ~ p/q, i.e. |A| — cq, 

= cp, with p and q positive integers relatively prime. 



8.1 Eigenvalues having the same sign 

We should first of all notice that for A/i > 0, no resonances are actually possible 
unless one of the eigenvalues is a multiple of the other, and in this case we have 
only one resonant term. 

Indeed for A/i > the only possible resonant terms are given by 

V = (^y^l^^^ (A//ieN) or V - (^J^^^ (/i/AeN). (8.2) 

In order to see this, recall the resonance relations are now miA + m2/i = A 
for the X component, and miA + m2/x = /i for the y component. These give 
(mi — 1)A + TO2/i — for the x component, and miA + (m2 — 1)// = for the y 
component; here mi,TO2 are non-negative integers, and mi + m2 > 2. As A, /x 
have the same sign, the only possibility in the x case is mi = 0, and A = m2/x 
with m2 > 2. Similarly, in the y case it must be TO2 = 0, and thus /i = TOiA 
with TOi > 2. 

We have thus shown that for A/x > the standard NF for the case S4 is 
given, with a, /? are arbitrary real constants, by 

X = Xx + ay^ , y — fiy for A//i = fc e N, fc > 2 
X = Xx , y = ^iy + l3x'' for /i/A = fc e N, fc > 2 (8.3) 
X = Xx , y = ^J■y otherways . 

In each of these cases the PRF is trivial, i.e. it coincides with the standard NF. 
We will thus give no further consideration to the case A/i > 0. 



22 



8.2 Eigenvalues with opposite signs 

Consider now the rational case with A/i < 0. Assume 

X = cq , 11= —cp (8.4) 

with p,q positive integers, relatively prime (no common factor), and c ^ a 
real number (notice wc could have p = q = 1, corresponding to /x = —A). For 
the sake of our discussion, we could as well take c = 1. 

The resonance relations give now {m\ — 1)A = — m2/x for the x component, 
and niiX = (1 — m2)/i for the y component. Hence we must have, for the x 
component, 1112/ {mi — 1) = — A//x = q/p, i.e. 

mi = kp+1 , m2 = kq . (8.5) 

Similarly, for the y component we have (m2 — l)/mi = — A//i = q/p, and 
therefore 

mi = kp , m2 = kq+1 . (8.6) 
Thus, the resonant vectors are of two types: 

Vfe - J ' - [{xPy'')''y) ■ 

Correspondingly we consider vector fields 

$fc = [{xPy''fx]d, , *fc = [{xPy'')''y]dy . (8.8) 

The most general NF will be in the form W = Wq + X](ci^^*fe +cif^*fe)) i-e- 

00 00 

X = Xx + J2 ci'\xPy'^fx , y = ny + Y^ c^k\^''y'fy ■ (8-9) 

The corresponding vector field will be denoted as W; its linear part is given by 

Wo = g$o -p*o- 

For our discussion it will actually be more convenient to consider linear 
combinations of the ^f/;, defined as 

= (2^) ((Z^fe-P^fc) = 2^5 - raa^) ; 

with this notation, the linear part Wq of the vector field W corresponds to 
Wo = 2cpqYo := (Yq. 

We also rewrite the corresponding vector field W, in view of the use of the 
vector fields Xk and Yk and for further discussion, as 

00 

W = CYo + J2(akXk + bkYk) (8.11) 
fc=i 
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where 

a, = (pc«+?cf) , h - (pcW-gcf). (8.12) 

Remark 12. Notice that $fe,^'ft G Wfe, but, with 2; = p + g, we have 
instead Xk,Yk e Wfc^. © 

The vector fields $fe and ^'^ satisfy the commutation relations 

[$fe,$m] = 7j(m-fc)$fc+m 

[4'fc,^'m] = q{m-k)'^k+m (8.13) 
and from these it follows that 

[Xk,Xm] = {m-k)Xk+m 

[Xk,Ym] = mYk+m (8.14) 

[Yk,Yrr,] = 0. 

Notice that these are the same as those encountered in discussing the case S3: 
we have thus to deal again with the Lie algebra Q = X y. Thus, provided 
we take into account remark 12, the algebraic computations considered there 
will immediately apply to this case as well. This correspondence between cases 
S3 and S4 is, of course, the one discussed in section 3 above. 



9 The S4 case: Poincare renormalized forms 

As remarked above, the algebra Kcr(£o) is spanned by vector fields {Xk,Yk} 
(with /c e N) which generate the same Lie algebra Q = X y encountered 
in discussing the case S3, as also discussed in section 3. The fact that the 
linear part is now given by CYq = '2'CpqYo, rather than simply by Yo, has no 
consequence on the discussion of linear subspaces, since the constants c, p, q are 
all nonzero and thus C 0: see remark 10. 

We can thus just repeat the discussion conducted in case S3, modulo remark 
12 above; we write again z := p + q. 



9.1 General results 

We will thus consider the terms in Wz, given by 
and consider the different cases 



ai ^ , 61 = 
ai = , 61 ^ 
ai ^ , 61 ^ 
, Oi = , 61 = 



(a) 
(b) 
(c) 
(d) 



(9.1) 



(9.2) 
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Notice that, considering a system which is already in standard normal form, 
the operators £i, ...,C,z-\ vanish; the first nontrivial higher homological opera- 
tor is 

:=[iy„.] = [aiXi + 6iyi, .] • (9.3) 

With exactly the same argument as in the discussion of the case S3 we have the 
following results. 

In case (a), where W^, = aiXi, Ker(jVli), Ker(jV1z_i) do just coin- 
cide with Ker(£o) (that is, the whole space on which the trivial operators 
A4i, ...Mz-i are defined), while Ker(A^2) reduces to the linear span of Wq 
and Wi, i.e. of Yq and Xi. The range of A4z is the whole linear span of the 
X, y, except the subspace spanned by Xi, X2, Yi. As vector fields Zi, Zi which 
are one in Ran(A1j;) and one in Ker(A4i) commute, no further normalization is 
possible. 

Thus, we obtain the same result as in case S3(a), with the role of A^i now 
effectively played by (which is the operator associated to Xi, and more in 
general to a\X\ + 61Y1). 

In cases (b) , (c) we do similarly reproduce the discussion of the correspond- 
ing cases of S3, again with the role of M.\ now effectively played by M.z- 

Thus we obtain the following expressions for the PRF when Wz 7^ 0: 

= C^o + 01X1 + 02X2 (a) 

< W = C^o+&iyi+Er=2Sfc^fe (b) (9.4) 
W = Cio + (ai^i + 6in) + a2X2 (c) 

In case (b), the ^-adapted procedure described in section 3 and subsection 6.3 
will actually give a more reduced NF (the LRF), see below. 

In the degenerate case (d) we should again proceed as in case S3(d). With 
the same meaning for as there and the same splitting in subcases (da), 
(db), (dc), we would obtain exactly the same expressions for the PRF as in 
(6.12), (6.16) and (6.20), at the exception of the linear part being given by (Yq 
rather than by lo- 
in case (db), using the ^-adapted procedure one would get, see (6.22), 

M 

W = CYo + a^X^ + 02^X2^ + J2 ^kYk ; (9.4') 

k=i' 

this also applies to the case (b), with v =1. 
9.2 Explicit computations 

The results of the explicit computations performed in the case S3 would also 
extend to the present case. Indeed, once we have transformed the original system 
into standard NF, the linear part Wq does not enter in the PRF algorithm any 
more, and thus the presence of the constant C (instead than one) cannot affect 
the computations in any way. Again, when translating the results obtained in 
case S3 to the present case, one has to take into account remark 12. 
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We will thus just follow the first steps of the computation in case (a) to 
illustrate this. 

We start from a system W^^^ which has already been brought to standard 
NF. Transformations generated by ft-fe G .Ffe for < /c < z are necessar- 
ily trivial, as for such k we have Ker(£o) ri = {0}, i.e. the Lie-Poincare 
"transformations" reduce to the identity. With a transformation generated by 

(z) 

hz = aiXi + (iiYi the term Wz is unchanged, while WjV^ is taken into 

W!f+^^ - 02^2 + {b2-aiPi)Y2 (9.5) 

and we can of course eliminate the I2 component by choosing [3i = b2/ai] we 
also choose ai — 0. 

Transformations with hk, z < k < 2z are trivial; thus we have = Wm'^^^ 
for all TO > 0. With a transformation generated by h2z — 012X2 + /32I2 (so that 

{'2z') 

h2z & Ker(£o) r\J'2z) the terms Wk, k < 3z, are unchanged; the term W^^ is 
taken into 

Wt^'^ - [«3 - ai«2] X3 + [(aifcs - ^2^2 - 2alf32)/ai] Y3 . (9.6) 
This can be eliminated by choosing 

a2 = as/ai , /32 = (0163 - 0262)7(201) . (9.7) 

Again, the transformations with generator hk G Kei{Co)f^J'k are necessarily 
trivial for 2z < k < 3z, and thus = for all to > 0. 

The explicit formulas obtained can be compared with those of section 6; they 
show the complete correspondence with the subcase S3(a). We believe there is 
no need to give further explicit formulas for the present case S4, as they can be 
read from the corresponding ones for case S3. 

Notice that also the expression of the PRF in terms of the vector fields 
Xk,Yk will be (except of course for the linear term Wq, where the constant C 
appears) exactly the same as in the case S3. 

10 Summary of results for other cases 

In this section we briefly recall, for the sake of completeness, the results obtained 
in ||3^, |3^ for the other cases where the PRF is nontrivial. These are cases S2 
and N2 of our basic classification of linear parts. An error contained in |^, ^ 
for one degenerate 82 subcase is also corrected. 

10.1 PRFs and LRFs for the case S2 

In case S2 one could work in C^, but we will stay within the framework of the 
present discussion and work in R^, i.e. we will deal with real matrices, and thus 
write 




(10.1) 
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Notice that in order to map this case into S3, it suffices to pass to polar coor- 
dinates. It is well known that, with — ~\- , the standard NF is then 



,2..^..^„.„^ (10-2) 



If the system is hamiltonian, then all the bk are zero; conversely, if all the 
bk are zero, the NF (2) is hamiltonian. Further reduction of the NF in this case 
has been studied by Siegel and Moser Q a long time ago (their results have 
recently been shown to generalize to higher dimensions |2§|] ). We will consider 
the general (non-hamiltonian) case. 

Let fi > 1 he the smallest k such that Ofe ^ 0, and ly > I the smallest k such 
that bk ^ (we assume both fj. and are finite). 

If /i < then (see subcases (a) and (da) for S3) the FRF is given by 

y — X + a^r'^^y + ar^^y . 
li ji — then (see subcases (c) and (do) for S3) the PRF is given by 

i = -2/ + r'^^'ia^x - b^y) + ar^^'x (10 A) 

y — X + r^'^{bf^x + a^y) + ar'^'^y . ^ ' ' 

Here ^ and 6^ ^ are the same as in (2), and the coefhcient a is a real 
number. A detailed proof of this result is contained in section 12 of a 
shorter proof is also given in Q . 

If !/ < /X, then (see cases (c) and (do) for S3) the PRF is given by 



X = -{l + b,r^-)y + EZf.r^'^akX 



In this same case, the LRF is given by 



X = -y + a^r'^'^x + ar^f^x - Yl,k=u V 
y = X + a^r'^'^y + ar^f^y + Y.k=y^kX. 



(10.5) 



(10.6) 



Here 7^ and bk (for v < k < ji) are the same as in (2), and the coefficient 
a is a real number. 

The computation for this case given in |3l| contained a mistake: the 
coefficients bk cannot be changed (via a PRF-likc transformation) to eliminate 
the corresponding rotation terms without producing radial terms. It should also 
be stressed that the reduced NF obtained in these papers, even after correction of 
this mistake, is the LRF (and is not obtained with the generic PRF procedure); 
in particular, in the case v < fi this is not a PRF according to our definition. 

It should be noticed that these results can be obtained, i.e. the case S2 
can be studied, more easily using the approach of the present paper, as we now 
briefly indicate. 
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We define, as in g^, dilation and rotation linear vector fields 

D = xdx +ydy , i? = -ydx + xdy (10-7) 

and with this compact notation, writing also (x^ + y^), we define 

^-fc := r^'' D , $fc := r^'^ i? . (10.8) 

It is immediate to check that these vector fields satisfy the commutation 
relations 

= 2(m - fc)*fc+™ , [$fc, $„] = , $m] = 2TO$fc+„ . (10.9) 

That is, we have the same algebraic structure as the one encountered in ana- 
lyzing previous cases. We can make it identical, including coefficients, see (5.3), 
by defining 

Xk = (1/2)1/3 , Ffc = (1/2)1/3 $fc . (10.10) 

In this way, the explicit computations performed for the case S3 can immediately 
be applied to this case as well. We write now the standard NF as 

oo 

W = CYo + ^(afeXfe + fefe^) (10.11) 

fc=0 

where C — (2)^/3. Formulas for the FRF can be read off the discussion of the 
S3 case. In particular, (6) corresponds to (6.22); see also (7.33). 

10.2 PRFs for the case N2 

In the (nonregular) case N2, we have 

and the standard NF is given by 

y = EZA^^x'^'-' + bk^'y) ^ ' 

Let ^ > 1 be the smallest k for which {aj. + bf.) 0; then the PRF is given 

by 

i - y + b,x^+' + ax^+' + Y.Z,+2 bkx'^' no 14^ 

y= a^x^^+^ +b^xf'y + axf'+'y + Y.T=ii^kxhj + bkx''+^) ^ ' ' 

This represents a very poor simplification of the standard NF; such a poor 
performance of the algorithm is related to the vanishing of the semisimple part 
of A, i.e. to the fact the singular point is nonregular. A detailed proof of (14) is 
contained in section 15 of |3^. For a discussion of this singularity, see |5^ and 
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11 Conclusions 



We have studied vector fields in the plane around a singular point by means of 
normal forms theory, discussing in detail all possible cases when the singular 
point is a regular one. In doing this we have assumed the linearization of the 
vector field has been preliminarly taken into Jordan normal form. 

We have also shown that when the standard normal form (NF) is nontriv- 
ial, the Poincare renormalized forms (PRFs) approach permits to substantially 
simplify the expression of the vector field in normal form. 

Thanks to constraints on the structure of the infinite dimensional Lie algebra 
of two-dimensional vector fields in normal forms p7| , |53| , there is a substantial 
correspondence between different cases where the NF is nontrivial, and compu- 
tations performed in one case can be mapped into any other one. We have taken 
advantage of this property, and performed explicit and detailed computations 
in one case (S3), using them for other cases as well. 

Considering the Lie algebra structure of vector fields in normal form also 
allow to define a different reduction scheme, designed to take advantage of this 
structure. The reduced normal form thus obtained, and called Lie renormalized 
form (LRF), is not necessarily a FRF; actually we have seen that in some of 
our subcases ~ i.e. cases (b) and (db) for all the semisimple linear parts - the 
LRF is finite while the PRF is infinite, and the LRF is not a special instance of 
PRF. The LRF approach is directly related to Broer's approach ^ and to 
Baider's work [Q. 

The local behaviour of vector fields in around regular singular points 
is of course very well studied, so that the real interest of our discussion is not 
in the expressions obtained themselves. Rather, we have shown that the PRF 
approach is viable to obtain a very explicit description of further reduced normal 
forms, even with the use of limited computing facilities. 

The detailed analysis given here also led to implement considerations based 
on the Lie algebraic structure of vector fields in normal forms, and to define a 
5-adapted reduction procedure, the LRF procedure, conjugating Lie-algebraic 
considerations a la Broer and the PRF algorithmic approach. This is of much 
wider use than the limited one considered here. 

We have also corrected a computational error contained in previous work, 
and clarified (see also the appendix) some confusion on the issue of PRFs present 
in the literature. 
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Appendix. The Bruno system 



As mentioned in remark 2 at the end of the Introduction, in his reviews of my 
papers |3^, for Mathematical Reviews , and again in his recent book p6| 
(section V.22) and in a preprint which he was so kind to send me, A.D. 
Bruno has claimed that the main resuh of my works js^, ^ is wrong; he also 
gave a "counterexample" to my result. 

This claim seems to originate in a misunderstanding of my previous papers, 
and the purpose of this appendix is to clarify the situation, also by completely 
explicit computations, dealing precisely with the system proposed by Bruno. 

The "counterexample" mentioned above falls in subcase S3(b) of the clas- 
sification considered here and was given in 113, U& as 



y = y + xy + x^y = {1 + x + x'^)y ; 
according to Bruno , the PRF for this system would be given by 



X = x^ 

y = y + xy = {l + x)y 

with no higher order terms. In this is changed to 

X — a2X^ + ax^ 

y = y + f3xy = (1 + /3x)y , 



(Al) 



{A.2') 



{A.2") 



again with no higher order terms, where a2, a, f3 are some constants (no mention 
is given to this discrepancy in ; in both cases no computation is reported to 
explain how these are obtained). More recently, Bruno and Petrovich also 
considered a slightly generalized form of this "counterexample" , see below. 

The discussion of section 5, and the explicit computations of section 6, show 
that (A.2) - in either one of its versions - is not the PRF for system (A.l). 

Actually, in his reviews, book, and preprint, Bruno quotes my result in a 
form which does not correspond to - and is not equivalent to - the statements 
I gave in |3l); thus his sweeping assertion that "it is easy to see that the 
statement of Gaeta is wrong" (see p. 275, [|l^) does refer to an incorrectly 
reported versio ngof my result. 

The key difference is given by the fact that the role played by the Aik 
operators in my construction is, in the version reported by Bruno ^6|, p^ , 
taken by the Ck ones (with no restriction to kernels of Cs with s < fc; see 
section 2). Thus, the result "quoted" by Bruno should not be attributed to my 
papers (incidentally, the statement attributed to me in [l5| , 16, is wrong in 
an obvious way). 

To avoid any confusion about PRF for the system (A.l), let us perform the 
Poincare renormalization algorithm up to terms in Wg. I stick to the proper 

^Notice that according to Bruno's definition (but with the notations of the present paper) 
the PRF would be characterized by the property that C'^_(Wk) = for j < fc in the first 
review , and for j = fc — 1 in the second of jis) and in ]l6| . 
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general PRF procedure as stated in |3l|, i.e. not consider the LRF scheme 
(actually the system (A.l) is already in LRF). I will freely use the notation 
introduced in discussing the case S3. 

We write W = x^dx + y(l + x + x^)dy. The system is already in NF, so 
we write W^^'' = W; the system is taken into PRF by operating successive 
transformations with generators Hk = akXk + PkYk (to avoid any possible 
misunderstandings, let us specify there is no sum on k). 

It results that choosing /3fc — and with 

ai = -1 , (12 = 1 , as = -2 , = 9/2 , , , 

as = -12 , ag = 33 , ar = -99 ^ ' 

the system takes successively the forms (where 0(9) denotes terms in Wg and 



higher) 














Yo+Yi 


+ X2- 


-X3-Y3+Xi + 2Yi+ 








-X5 -3Y5+X 


e+AYe-Xj- bYj + Xg 


+ 6^8 + 


0(9) ; 


M?(3) = 


Yo + Y, 


+ X2- 


-X3 + Xi + 2Yi - 2X5 - 


(9/2)n- 


f 




+3^6 + 




(11/2)X7 -25^7 + 9^8- 


h 5AYs + 


0(9) ; 




Yo + Yi 


+ X2- 


-X3 + Xi- (9/2)n+ 








+3^6 + 


UYf,- 


(15/2)X7- 33^7 + 13^8 


+ 99Ys ' 


fO(9) ; 


M?(5) = 


Yo + Yi 


+ X2- 


-Xs + Xi- 6Xe+ 








+l2Ye - 


3X7 - 


331^7 + 13^8 + 99y8 + 0(9) ; 






Yo + Y, 


+ X2' 


-X3 + Xi- 6Xe + 33X7- 


f 






-33Y7 - 


nxs' 


K 99y8 + 0(9) ; 






= 


Yo + Y, 


+ X2- 


-X3 + Xi- 6X6 + 33X7 


- 143^8 


+ 99Y8 




Yo + Y, 


+ X2- 


-Xs + Xi- GXe + 33X7 


- 143^8 


+ 0(9) 



0(9) ; 
(A.4) 

The latter is the PRF, up to terms 0(9), for (A.l). 

In ||l7|| , Bruno and Petrovich consider general systems with linear part cor- 
responding to our case S3. The standard normal form for these is, as discussed 
above, of the form 

k=l ( A 

y = y + Efe=i hx'^y ; 

we denote by m > 1 the smallest k such that Ofc 7^ 0, and by ^ > 1 the smallest 
bk such that bk 0. That is, we assume = for k < m, and bk — for k < i. 

Bruno and Petrovich suggest to consider the case £ < m < 00, and give as 
an example the case ^ = 1, to = 2. They claim, see formula (5.2) of ||l^, that 
the PRF in this case is 

y = y + px^y ; 
with Qm 7^ 0, and a, (3 real coefHcients. 
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Once again this formula does not correspond to the PRF computed in above, 
so that their proof that the original system (A. 5) cannot be conjugated - even 
formally - to system (A. 6), does not concern PRFs. 



Let us now consider the domain of convergence of the normalizing transfor- 
mation explicitly considered above. Here /3fc = for all /3, so that the mappings 
do not act on y. It is easy to obtain (with the help of an algebraic manipulation 
program) explicit expressions for the changes of coordinates and thus for the 
domain of analyticity of the overall transformation up to step k. 

We write the combined effect of the first k changes of coordinates as x ^ 
7j.(k) _ x/Bk{x); the denominators Bk{x) can be written in recursive terms as 

Bk{x) - [{Bk-i{x)f - (-l)Sfca;']'^'' • {A.7) 

The first numbers of the sequence 7^ (fc — 1, 2, ...) are given by 1, 2, 6, 18, 60, 
198, 693; obviously Bq{x) = 1. We will omit the derivation of this recursion 
formula. 

This also allow to explicitely determine the domain of analiticity of the 
transformations, which can be read from the roots of -Bfc(x); but the expression 
so obtained quickly become extremely involved and of little interest. 

(k) (k) 

Thus I have computed analytically the x_ ,x\ such that the overall trans- 
formation up to step k (i.e. x x^'^^) is analytic in the strip x^^* < a; < x^J^\ 
but just report here their numerical value. These are: 

x'l^ = -1 , x'l^ = -0.333333 , x'-^^ = -0.270929 , x''^^ = -0.244594 , 
x'-^^ = -0.228796 = -0.21915 , xL^^ = -0.21224 ; 
a;^^^ = 00 , a;^^^ = 1. , a;^^^ = 1. , a;^^^ = 0.668534 , 



x*/^^ ^ 0.668534 , a;^®'o.561419 , a;^^^ = 0.561419 



(A8) 



Finally, let us clarify the relation of the computations presented in this 
appendix with the problems, mentioned above, concerning the discussion of 
case S2(db)) with < ^ in ||^, |l|. 

It should be remarked that if one maps (as discussed in section 3, or simply 
passing to polar coordinates) case S2 to case S3, then the wrong result given in 
[ pO| |3l[ for the degenerate case (with the notation of the present paper, S2(db)) 
with z/ < /i would map exactly to the "pseudo-PRF" (A. 2') given by Bruno in 
[|5j . Notice however that there it is claimed that the error lies with the general 
statement and not with the computations of the example; also, (A. 2') is not 
claimed to be derived from the results for S2 given in [^o|, |3l[, but just to be 
the PRF according to the definition reported there. Thus, such a statement 
appears somehow mysterious; however, as discussed above, it does not involve 
PRFs according to their definition given in js^, ^ and in this paper, so that 
we don't have to deal with it. 
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